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1. Introduction 


HEXAMETHYLENE-TETRAMINF is one of the few organic substances crystallising 
in the cubic system. Though its molecules contain as many as 22 atoms, 
they possess tetrahedral symmetry. The study of the Raman spectrum 
is therefore of exceptional interest. The substance has been previously 
investigated by Krishnamurthi (1931); and Kahovec, Kohlrausch, Reitz and 
Wagner (1938). The polarisation of the Raman lines has, however, not 
been determined so far. In the present investigation, an intense photograph 
of the spectrum of a saturated aqueous solution of the substance has been 
obtained with a spectrograph of high dispersion and high resolving power. 
Some hitherto unobserved detail has been recorded and the polarisation 
character of the lines has also been ascertained. The results are discussed 
in the paper in relation to the structure of the molecule. 


2. Experimental Arrangements 


The substance used in the present investigation was Kahlbaum’s pure 
chemical and was fairly free from dust. A saturated solution of this 
substance in double distilled water was rendered dust free by carefully 
filtering it several times and was finally transferred into a double-walled 
Wood’s tube and sealed off. A few drops of concentrated nitric acid was 
added to the solution to render it free from fluorescence. The outer 
iacket of the Wood’s tube was filled with a concentrated solution of 
sodium nitrite solution which reduced the strong continuous background 
accompanying the Raman spectrum by cutting off all the ultra-violet 
radiations upto 4046 A.U. Two pyrex mercury arcs placed on either side 
of the Wood’s tube and in close contact with the latter served as intense 
sources of illumination. A Littrow glass spectrograph having a dispersion 
of about 10 A°'mm. in the 4358 A° region, and a resolving power of about 
3cm.- in the region between 4046 A° and 4358 A° was employed to obtain 
the Raman spectrum while the polarisation picture was taken with a Hilger 
two-prism glass spectrograph of high light-gathering power. 
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3. Results 


The Raman frequencies of hexamethylene-tetramine as well as their 
polarisation characters obtained in the present investigation are listed in 
Table I along with the results reported by Krishnamurti and Kahovec 
et al. 











TABLE I 
Raman Frequencies of Hexamethylene-Tetramine 
No | Author | Polarisation Krishnamurthi Kahovec et al. 
| | 
1 265 (4) (K.e) § D 260 (1) (e) 
2 299(1) ,, | D 284 (1) (+42) 
3 398 (4), D od bis 
4 455 (4) ,, | D 457 (m) 454 (4) +(e, c) 
5 514 (2) D 513 (S) 508 (3) +(e, c) 
6 68444) ,, | D sh 682 (4) (e) 
7 782(10) ,, | p 782 (st) 782 (10) (f, e, ¢) 
8 821 (4) | D 820 (S) 818 (0) (e) 
9 1009 (3) ) D 1008 (m) 1011 (6b) (e) 
10 1023 (1)) D ene ie 
11 1048 (6) P 1052(m) 1046 (6b) (e) 
12 1093 (1) D ne 1090 (2) (e) 
13 1240 (2) D at 1243 (2) (e) 
14 1346 (4) D 1346 (m) 1345 (5) (e) 
15 1384 (2) D ce 1377 (2) (e) 
16 1440 (1) ) D 1452 (m) 1451 (5b) (c) 
17 1452 (6)} D ai ae 
18 2890 (2) P mi 2873 (3) (e) 
19 2965 (2) v 2948 (st) 2961 (1) (e) 
20 2973 (2)) D wt a 
21 2989 (2)) D 2977 (st) 2980 (3) (e) 











One new line at Av = 398 has been obtained in the present investiga- 
tion in addition to all the frequencies reported by Kahovec et al.; while 
some of the lines reported by these authors as broad have been found to be 
doublets, e.g., 1011 (65), 1451 (55) and 2980 (3). On the other hand the 
line at 1046, reported as broad by the above authors is found to be quite 
sharp and intense. These facts can be clearly seen from the spectra 
reproduced in Plate II. 

4. Discussion 

Gonell and Mark (1923), Dickinson and Raymond (1923), and Wycoff 
and Corey (1934) have investigated the crystal structure of hexamethylene- 
tetramine. The unit cell is found to contain two molecules of C,H,.N,. 
The molecule C,H,;N, possesses a symmetry Td with atomic distances 
C-C = 2-58 A° and C-N= 1-48 A°. The carbon atoms thereby occupy the 
corners of a regular octahedron as seen in Fig. 1. The N atoms which 
occupy the corners of a tetrahedron lie on the threefold axes of symmetry 
of the octahedron. The character table, the number of modes of vibration, 
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Fic. 1. Structure of Hexamethylene-tetramine molecule C,H,.N, 


and the selection rules for the Raman, and Infra-red frequencies are given 





















































in Table II. 
TABLE II 
Character Table for C,H,.N, 
| Total Internal Trans- | Selection 
number | vibrations lations rule 
Td E |&8C,|3C,| 60 | 6S, of and 
. vibra- ; Rota- 
is tions | Mole-/| Ring | tions | Raman | Infra- 
| cule Red 
e | = 
a4 A, Pelee eees o 4 4 a S P f 
e A, af 2] 0 7 0 la [-2 1 1 0 | Oo f f 
b «A tee 1 eieie 5 5 a oe D f 
€ F, | 31 O f—t |=! 1 7 6 2 | f f 
F, 1 3{ 0 j-1 | 1 [-1 | 10 9 4 IT D 4 
U, wy: IR SES SES 
hyd, ..|66| 0 |-6 | 36 | 0 | 
hip,’ ..(60| 0 | 0 |36 | 0 | | 
| | 
ff : , : 
5. Symmetry Co-ordinates for Ring Frequencies 
3 The principal axes X, Y, Z are obtained by joining the carbon atoms 
5 (5, 6), (7, 8) and (9, 10) the positive directions being also the same. The 
. Co-ordinates of any mode of oscillation are expressed with respect to a 
: parallel co-ordinate system with its equilibrium position as origin. The 
y following “normal co-ordinates are obtained :— 
, Ql yt Xe Xg— Ket Yi— Yat Va— Vat Am Za Za Za | 2A 
’ eal | 
Q X5— Xet Ya— Vet 29— 219 
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aa Vi-— Vat Va— Va— 24 + Ze t+ Zy— 2% 
Q'35 2X, + 2Xg— 2Xg— 2Xq— Vit Va— Vat Va— 2+ Ze Z3— 2, 
ta Yr—- Va- 29+ 210 
Qu 2X5— 2X—e— Yr + Ve— Zot Z10 
Que My? (i —Ya— Va HV +21 —Z2 +23 —24) +M ct (27-23 — Yn +Yy0) 
oo My? (x) —X2—X3 +44 +25 +22 —Z3 —Z4) HM ct (Xp —X 19 —Z5 +2¢) | Fy (1) 
Qin: Myt (4) —X2+X%3—X 4 t+Vi + ¥2—Ya—Va) +Mc# (x,—X5 +e —Ys) 
QVse 4 Mc'i—Y2— Ys +¥g +21 —22 +23 — 24) —6 My? (2, -23 — 9 + V0) 
‘36 4 Mch(xXy — X29 —X3 +4 +21 +Z2—2Z3—24) —6 My? (x99 —25 +2¢ F, (2) 
se 4 Mc H(xy —X2 +X3—Xa +V1 +¥2—Vs— Ys) —6 My? (x; —Xg +¥6—Vs 
Q. My (4 + X2+ X3 + X4) + Mc? (5+ Xe + Xz + Xg + Xy + X10) | 
Q, My? (¥i+ Yat Yat Yad + Met (¥5+ Yet ¥2+ Vet Vot Vo) | F, (1) 
Q Myt (4+ Zot 25+ 2) + Met (25+ Zot 27+ Zt Zot 20) | 
x 4 My? (x5+ X6 + X7+ Xg+ Xy4# X40) — 6 Met (x1 + Xg+ Xy+ x4) 
Q'S 4 My? (5+ Vet Ve + Yet Yot Vio)— 6 Mct (Vit Vo+ Yat Va) 
‘oc A My? (25-4 26+ 27+ Ze+ Zp+ Z10) — 6 Mct (21+ Zet+ 23+ 24) ! 
QVx2 Yi Va— Vat Va 241 Ze— Zyt 2 
Q'y Xam Xa Xa t Xq— 21— Ze t Zzat 2%, 
QV. Xam Ket X3— Xg— Vi— Vat Vat Va 
Q's: 27— 28+ Yo— Vo 
Q's  X9— Xo + 25— 26 
Q's X7— Xge— Ve + Vs 
The ten modes of oscillation of the ring are diagrammatically represented 
in Fig. 2. 
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Analysis of the Raman Spectrum of Hexamethylene-tetramine 


According to group theoretical analysis the sixty vibrational forms of 
the molecule fall into five classes given in Table II. Most of these vibrations 
become degenerate due to the high order of symmetry possessed by the 
molecule. Taking into account the degeneracy as well as the selection rules, 
18 frequencies are expected in the Raman Effect and 4 in Infra-red absorp- 
tion. Out of the 18 frequencies active in the Raman effect 8 are eharacter- 
istic of the ring and 10 are characteristic of the (C-H) vibrations. All the 
8 ring frequencies as well as the 6 deformational (C~H) vibrations ought 
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The Raman Spectrum of Hexamethylene-tetramine 
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to fall in the region below 1500 cm.-! while the remaining 4 (C-H) frequencies 
should appear in the region near 3000 cm.-? 


According to theory two of the ring frequencies ought to be polarised. 
These can be easily identified as 782 and 1048; they being the only two 
polarised lines below 1500cm.-' The other 6 depolarised ring frequencies 
can be identified with the lines 299, 398, 455, 514, 1009 and 1023 while 
the 6 (C-H) vibrations falling in this region can be taken to be 1093, 1240 
1346, 1384, 1440 and 1452. The four (C-H) frequencies expected in the 
region near 3000 cm.-! are 2890, 2965, 2973 and 2989 of which the former 
two are polarised and the latter two depolarised; which is in agreement 
with theoretical expectation. 


All the expected lines are thus accounted for barring the three weak 
lines at 265, 684 and 821. They may, however, be considered as the 
differential or combinational tones, as is the case, of the following frequen- 
cies :— 


265 = 1048— 782 
684 = 1452— 782 
821m 514+ 299. 
The observed Raman frequencies are thus found to be in good agreement 


with the tetrahedral structure of the molecule. 


In conclusion the author thanks Prof. Sir C. V. Raman for his kind 
interest in the work. 


Summary 


The Raman spectrum of hexamethylene-tetramine in saturated aque- 
ous solution has been obtained along with the polarisation character of the 
lines; the latter being obtained for the first time. One new line at 398 has 
been observed in addition to all the frequencies reported by Kahovec et al. 
The lines at 1011 (6 5), 1451 (5 6) and 2980 (3) reported by the above authors 
as broad have been found to be doublets; while the line at 1046 reported 
as broad is found to be quite sharp and intense. The number of distinct 
frequencies recorded and their polarisation characters show satisfactory 
agreement with the tetrahedral model assigned for the molecule from X-ray 
studies. 
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More recently Sir C. V. Raman (1941, 1943) has reconsidered the 
problem of the lattice vibrations and has shown that the experimental facts 
of spectroscopy necessitate a completely new approach to this problem, 
His theory has been successfully applied by the present author to the 
evaluation of the specific heats of rocksalt, diamond, and the metals 
(Dayal, 1944 a,b,c). In the present paper the calculations have been 
extended to KCI for which an extensive specific heat data is available, and 
to KBr, the only other crystal of this class which possesses such data. 


2. The Frequency Spectrum in the Born and the Raman Theories 


Let us for simplicity take the case of a general cubic lattice, the unit 
cell of which is described by the three equal orthogonal vectors a,, a, a, 
The cells are distinguished by the cell index / where / stands for the three 
integers /,, /, and /,. The p different atoms of the cell are indicated by the 
particle index k. The equilibrium position of the Ath atom in the /th cell 
is given by r(k, 1) = la, + /oa,+ /a,+ r (k, 0) where r (k, 0) is the equili- 
brium position of the particle in the cell at the origin. According to Born 
(Born and Goppart Mayer, 1933) the displacement ui, of the (/, k)th particle 
is given by the plane waves of the type 
w =U; it 2avt + [Sr (k,D]} (I) 
The direction of the vector S is the direction of propagation of the wave 
and its magnitude is equal to 4 where A is the wave-length. The possible 


values of S are determined by the cyclic postulate, according to which the 
infinite lattice is subdivided into cubic blocks of N® cells, N being supposed 
to be a very large number. It is assumed that the vibration pattern repeats 
itself exactly in each of these blocks so that the vector S is given by 


= ((% {le {%s ) 
where the b’s are the lattice vectors inverse to the a’s and A, Mo, ng are 


integers given by the condition — . <n< “4 If the scalar components 


of the wave-number vector 1 have the values I : a ; ; 
A A, As’ Ag 


that A; = “24 where a is the lattice spacing. For the waves travelling along 
the cube axes A; are directly equal to the phase wave-lengths and are the 
sub-multiples of the length of a side of the large cube of N? cells. 
Because of this assumption, they bear no relation to the crystal spacing and 
crowd themselves increasingly as they became smaller: It can be easily 


it follows 
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sen that one half of them (4th from each branch, optical and acoustic) have 
values intermediate between 2a and 4a, } between 4a and 6a and so on, the 
proportion found in any given range of wave-lengths diminishing continu- 
ously and vanishing ultimately for great phase wave-lengths. 


The nature of the frequency spectrum in the theory of Born is essentially 
determined by the cyclic postulate. There are two sections of it namely, 
the acoustic and the optical branches. The frequencies in the acoustic 
branch increase continuously from zero which is the value for infinite phase 
wave-length, to a limiting value in the short wave-region. In the optical 
branch, the frequency for great phase wave-lengths has a value approximately 
independent of such wave-length but generally falls off more and more 
rapidly with the diminishing wave-length, till it reaches its minimum value 
at the short wave limit. The optical frequencies thus form a continuous 
spectrum extending between the limiting values of the frequencies corres- 
ponding respectively to the infinite and the smallest permissible wave- 
lengths. [Illustrations of its nature will be found in the drawings for the 
rocksalt crystal given by Kellermann (1940). 


On the other hand, the theory of Raman starts from the consideration 
of the theoretically possible normal modes of vibration of the atoms about 
their positions of equilibrium in the structure of the crystal. Alike in the 
case of polyatomic molecules and of continuous systems, the fundamental 
property of a normal vibration is that all the particles in the system have 
at any instant the same or the opposite phase of vibration (Rayleigh, 1926). 
Hence, each of the normal modes with which we are concerned must satisfy 
this condition and also be consistent with the fundamental property of the 
crystal, namely, its translational symmetry along each of the three axes of 
the Bravais lattice. Raman has shown that there are (24p—3) modes of 
atomic vibration in a crystal which satisfy these requirements. The atomic 
vibration spectrum thus consists of (24p—3) discrete or monochromatic 
frequencies, the number, however, being further diminished by degeneracy 
in the crystals possessing types of symmetry other than that of the trans- 
lational kind. In these (24p—-3) modes, equivalent atoms in the crystal 
structure possess the same energy of vibration, but the phase of the vibra- 
tion is either the same for all such atoms along an axis of the Bravais lattice 
or else is alternately positive or negative. It follows that the (24p—3) modes 
represent the possible internal vibrations of the atoms contained in a 
space-unit having 8 times the volume of the lattice cell of the crystal, each 
of the modes having a degeneracy '/,th as large as the total number of lattice 
cells in the crystal. This degeneracy may be regarded as a consequence of 


the transaltional symmetry of the crystal structure. 
A2 
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Chelam (1943) has discussed the modes of vibration of the NaCl type 
of lattice on the basis of the Raman dynamics. His results are shown in 
Table I. It will appear from this that the possible different modes are just 
those which in Born’s theory appear as the limiting vibrations with the 
greatest and the least possible phase wave-lengths normal to the faces of 
the cube and the octahedron. Lyddane and Herzfeld (1938) have described 
the various frequencies of the Born theory by a vector + the scalar compo- 


9 
nents of which parallel to the cube axes are given by 7; = i For conveni- 
z 


ence of reference we have given the values of 7; corresponding to the Raman 
frequencies below. Kellermann has used the quantities p,, p,, p, in this 
connection which can be obtained from the 7; on multiplication by twenty. 


TABLE I. Modes of Vibration of the Lattice Structure 





Designation and Description of the Mode Degeneracy | rs, Ty, T, 





vy, Oscillation of the alkali and the meee lattices against 
each other ce ‘ ot 3 0, 0, 0 


ve Oscillation of the alkali layers against each other normal 
to the octahedral planes with the halogen layers at rest 4 4, 4.4 


vg Oscillation of the halogen layers against each other 
normally to the octahedral ae with the alkali atoms 





at rest .. de : x es if 4 4, 3.4 
vz Same as vz but transversely oe st me ry 4, 4,4 
vs Same as v3 but transversely ese ou Pe 8 4, 4, 3 

vg and y, Coupled oscillation of the alkali and halogen ions normal 
to the cube planes ay: Bet = 3 4,0,0 

vg and vy Coupled oscillation of the alkali and the a ions 
transverse to the cube planes Fie : . 6, 6 4, 0,0 
Translations .. * bss “4 3 000 

Total ‘sa 48 











3. The Evaluation of the Frequencies 


Since the discrete frequencies given by the Raman dynamics are included 
also in the continuous spectrum of the Born theory, it is possible to use the 
general methods devised by Kellermann (/.c.) and others for their evaluation. 
The displacement of an atom is again given by the equations (1) and (2) 
but 7 is restricted to the values shown in Table I against the Raman 
frequencies. We have generally followed Kellermann’s treatment in the 
present paper. It is supposed that the potential energy ¢ between any two 
atoms depends only on the nature of the atoms and their distance of 
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separation. For the atoms (/, k) and (/’, k’) it may be written as 4,3”. The 
displacement components of the atom (/, k) parallel to the axes are 
denoted by ui,, ui, and ui. respectively. If the oscillations are harmonic, 
the force due to the mutual displacement of any two atoms is given by the 
product of the appropriate force constant and the relative displacement. 
The force constant being a second differential coefficient of the potential 
energy, the X component of the force on the atom (/, k) due to the displace- 
ment of the atom (/’, k’) is given by 
a 


ad ahd , 
—— (i ~ al (3) 
y 


OX OVP 
where the summation is to be carried out for all the three co-ordinate axes. 
The total force on the atom (/, k) due to the displacements of all the atoms 
is written as 


m i, =2TZZ (Ge ey Uby (4) 
k’ qu’ y 
with the abbreviations 
‘ bh 
1-V i kh 
( kk Jey ve IXO¥ 


and 
(Pe ey = -22 (thy ey 
k’ 


when the dash over the summation sign indicates that the term /=0 and 
k=k’ (the action of the point on itself) is to be excluded. The deriva- 
tives are all to be evaluated for the equilibrium positions of the atoms. 
The displacement of each kind of particle is given by the equation (1) in 
which S; is to be obtained from the 7; appropriate to the different frequen- 
cies. On substitution of (1) in (4) we get 

47? v?m, Uz, + - & [kk': xy] Uy, +0 (5) 
where ee 

[kk’ : xy] = EZ) (Headey eS TH 6) 
r., being the vector joining the equilibrium positions of the (0, k) and the 
(I, k’) atoms. The coupling coefficients (6) occur as the coefficient of the 
amplitudes U;, in the equation (5). There are six such equations, and on 
eliminating the amplitudes from them we get a sixth order determinant for 
the frequencies. For all the modes of vibrations occurring in the Raman 
dynamics, the determinant splits up into factors so that the solution of the 
equation is very much simplified. There are mainly two kinds of forces 
acting between the atoms of the lattice. The first is the electrostatic force 
between the ions which, because of Coulomb’s law, decreases very slowly 
with the distance. The other is the short range repulsive force which keeps 
the lattice from collapsing, and is practically confined to the neighbouring 








82 Bisheshwar Dayal 


atoms. The summation in the definition of the coupling coefficients is carried 
out separately for each kind of force. That due to the electrostatic forces 
alone involves some difficulty and has been carried out for a number of 
frequencies both by Lyddane and Herzfeld (/.c.) and by Kellermann (i.c.). 


; ; ‘ e 
The numerical values have been given by them in the units of o where e is 


the charge on the ion, so that on substitution of the values of these con- 
stants the results can be applied to all the NaCl type of lattices. Except 
for the case *=0 the results of both the calculations are the same. For 
this vibration, Lyddane and Herzfeld have obtained different coefficients 
for the transverse and the longitudinal waves. Kellermann has questioned 
the validity of this result and has obtained only one value of the coefficients 
for all the three modes belonging to this case, the values coinciding with 
those given by Lyddane and Herzfeld for the transverse vibrations. We will 
follow Kellermann and use only one set of coefficients for += 0. 


The evaluation of the coupling coefficients due to the repulsive forces 
is not very difficult. Lyddane and Herzfeld have used the Born and J. Mayer 
(1932) expression for the repulsive potential and have thus omitted the small 
contribution of the Van-der-Waals forces towards the coefficients. Keller- 
mann has, however, taken this into account in a rough way by using the 
observed value of the compressibility and the equilibrium condition of the 
lattice. For the actual evaluation of the coefficients, Lyddane and Herzfeld 
have introduced two constants c, and c, given by 

— = 2 i r= a| 


4rre2 [dv (r) 
S tae Be ioe 
az © 2} or | “ a) 


(7) 


where v(r) is the repulsive potential of the ions. Kellermann’s constants 
A and B are 87 times c, and c, respectively. If we use Kellermann’s 
method, c, is the same for all the NaCl type of lattices and is equal to 
— +0466. c, has a different value for each salt and is given by 


__ 1 sla 4) 


where « is the compressibility at the absolute zero of the temperature, 
a the Madelung constant and e the ionic charge. From their calculations 
for the coupling coefficients, both the authors have obtained the following 
formule for the frequencies. For v,(7=0), only Kellermann’s expression 
is given. w stands for 47v and p is the reduced mass, 


(8) 
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Gre? 
w= as [cy + 2cg-- $] 


2 
wa? and wg? = “©, fe, + 2cq-+ +2876] (k= 1, oF 2) 
ke 
' (9) 


2 
ofutiat« oa [cy + 2¢q— 1438] (k= 1, or 2) 


aro (Ye 


where for w, and w, A= -1724+¢,+ 2c, 

B= -5984—c,+ 2c, 

and for w, and w, A=—-0862+ c,+ 2c, 
=— +2992+ c 


In the paper on the specific heats of rocksalt (/.c.), we pointed out that 
the frequencies obtained by Kellermann were slightly lower than those 
observed in the Raman effect. This is evidently due to the use of the room 
temperature compressibility in obtaining the constant c,. In the calcula- 
tions on KCI we have used the value of the compressibility at 80° K obtained 
from the measurements of Durand (/.c.). For KBr where the low tempera- 
ture measurements are not available. c, has been calculated directly from 
the experimental value of the infra-red frequency by the use of the formula 
for »,2 in the equations (9). For KCl we have used «= 5-03 « 10712, 


¢,="491 and “5 = 7-48 x10%. In the case of KBr c, = 506 and . = 6-4 








4 


x10%. The numerical values of the frequencies are given in Table II 
below. 


Tasie Il. Vibration Frequencies of KCl and KBr 





KCl KBr 





Frequency Designation Degeneracy 
| cm | Le Cm.} hv 
| | 
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It will be seen that the infra-red proper frequency », for KCl jis 
exactly the same as the one observed by Barnes (i/.c.). The use of the 
infra-red frequency for evaluating ¢, in the case of KBr seems, therefore, 
to be perfectly justified. 


4. Evaluation of the Specific Heats 


We have so far considered the nine optical frequencies which comprise 
15,16 of the total degrees of freedom of the lattice. In the earlier papers 
on rocksalt and diamond, we had assigned the 1/16 residual degrees of 
freedom of translation to the elastic spectrum and had evaluated their 
contribution to the thermal energy by means of a Debye function. The 
usual Debye characteristic temperature was, however, divided by ~/16 to 
allow for the reduced degrees of freedom. We will follow the same method 
here. The Debye characteristic temperature for KCl as evaluated by Durand 
from the elastic constants at 80° K. is 246, but the value which best explains 
the experimental data near 3° K. is 228. According to Blackman (1935), 
the two values should be the same. In view of the large uncertainty which 
is usual in the determination of the elastic data, we have used 228 as the 
correct Debye characteristic temperature. The reduced value is accordingly 
90. Blackman (1938) has given 179 as the value of the characteristic tempe- 
rature for KBr, from which we have obtained the reduced value as 71 for the 
purpose of these calculations. c,, are given by the formula 


C, (molar) = = [2 n,E (H) +3D (?)] (10) a 


1 





n, stands for the number of degrees of freedom pertaining to the frequency 
v,, While E and D stand for the Einstein and the Debye functions respec- — 
tively. The experimental data for KCl have been taken from Keesom and 
Clark for (3-17)° K. and from Southard and Nelson above that tempera- 
ture. The former experimenters have pointed out that on account of the 
small quantity of KCI used for the specific heat determinations, the accuracy 
of their results is not very high, while the latter have claimed an accuracy —_ 
of less than one per cent. C, have been taken from Venkateswaran’s paper 
(1941) who has given the expression 2c,*T for (c,— c,). For KBr, c, is taken Sir 
from Nernst and Lindemann (1911 /.c.). The last two readings are mean this 
values over an extended range of temperature. The calculated values of 
c, for KCl and KBr are given in Tables III and IV respectively. 
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TABLE IIT. Specific Heats of KCI in Cals. gm. mol. deg. 






































Temperature Optical | Elastic Total | Experimental 
°K frequencies | spectrum os 
| 
2-99 is -00215 -00215 | -00206 
3-52 om -00348 00348 -00338 
4-01 me -0052 -0052 -0042 
4-6 a -0077 -0077 -0072 
5-72 ne -0149 -0149 -0145 
6°52 Ry -022 022 -019 
7-65 8 -0353 0353 -0368 
9-2 -0028 -0601 0629 -0620 
10-06 -0062 -0767 -0829 -0840 
14-73 -087 -190 *277 +310 
17-09 186 +253 -439 -458 
21-21 -475 -349 *824 +842 
32-41 1-830 -521 2°351 2°358 
39-86 2-935 -585 3 +520 3-520 
49-27 4-278 -633 4-9] 4-91 
59-61 5 +548 *667 6°22 6-30 
69 -92 6°555 . 686 7°24 7+33 
79-11 7 +287 -698 7-99 8-12 
89-28 7-927 -708 8-64 8°77 
101-31 8-537 -716 9-25 9-35 
121-41 9-205 -724 9-93 10-04 
141-17 9-670 +729 10-40 10-60 
163-45 10-006 +733 10-74 10-96 
179-76 10-201 +735 10-94 11-13 
201 +32 10-395 -737 11-13 11-32 
240-7 10-63 74 | 11°37 11-58 
TABLE IV. Specific Heats of KBr in Cals. deg. gm. mol. 
| 
Temperature | Optical Elastic Total Experimental 
°K | frequencies spectrum G os 
| 
78-7 8-452 “715 9-17 9-42 
82-5 8-659 -718 9.38 9.46 
85-4 8-802 -719 9-52 9-58 
89-2 8-968 721 9-69 10-00 
137 10-133 734 10-87 10-66 
234 10-795 | -741 11-54 11-74 








In conclusion, the author expresses his grateful thanks to Professor 
Sir C. V. Raman, Kt., F.R.S., N.L., for his guidance in the course of 


this work. 





Summary 


The specific heats of KCl and KBr have been evaluated on the basis 
of the Raman theory of lattice vibrations. The nine discrete frequencies 
tequired by the theory and comprising 15/16th of the total degrees of 
freedom of the crystal have been calculated on the lines followed by Keller- 
mann and by Lyddane and Herzfeld. The electrostatic coupling coefficients 
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have been taken from the work of the above authors, while the repulsive 
coefficients have been calculated from the low-temperature value of the 
compressibility and the equilibrium condition of the lattice. In the case of 
KBr where the low temperature value of the compressibility is not available, 
one of the repulsive coefficients has been calculated from the known value 
of the infra-red frequency. The residual 1/16 of the degrees of freedom 
not covered by the optical frequencies has been allotted to the elastic 
spectrum whose contribution to the thermal energy has been evaluated by 
using a Debye function with a greatly reduced characteristic temperature. 
An excellent agreement between the theoretical and the experimental values 
of c, has been obtained throughout the entire range of temperature. 
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7. Introduction 


SILICON crystallises in the cubic system with a diamond-like structure but, 
contrary to the general belief about the applicability of the Debye function 
to the cubic crystals, its thermal energy is not capable of being explained on 
the Debye theory. Its specific heats have been measured over an extensive 
range of temperature between 20 and 1000 degrees absolute. The 
characteristic temperature which best explains the experimental data at high 
temperatures has been given by Magnus (1923) as 712, but as the temperature 
is lowered, the experimentally observed specific heat gradually becomes 
greater than the value calculated from the Debye function, and at 40° K., 
the former becomes four times as large as the latter. So far, the only attempt 
to explain this so-called anomaly has been made by Simon (1926) on the 
basis of Schottky’s quantum jump theory. He considers the specific heat 
as made up of two parts. The first is the normal part of the specific heat 
which is given by a Debye function (@ == 790), while the second is attributed 
to the presence of the two distinct atomic states in the lattice which are 
independent of the normal modes. The contribution of the latter to the 
atomic heat is assumed to be given by a Schottky function (U = 246). This 
theory has been criticised by Fowler (1935, 1936) who finds no a priori justi- 
fication for assuming the existence of two distinct atomic states in such cases. 
According to him, these specific heat anomalies arise from the fact that the 
Debye theory does not give the correct vibration spectrum of the solid. 


More recently Sir C. V. Raman (1941, 1943) has reconsidered the nature 
of the vibration spectrum of solids and has shown that it is completely differ- 
ent from the continuous spectrum given by the Debye theory. The present 
author has found (Dayal, 1944 a-d) that the specific heats of the alkali halides, 
diamond and the metals can be completely accounted for on the basis of 
the new theory. In the present paper it is proposed to extend the investiga- 
tion to the case of silicon. 


2. Evaluation of the Frequencies 
In the analogous case of diamond, the vibration spectrum is com- 
pletely known as a result of the spectroscopic investigations. All that is 
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necessary is to assign the frequencies to their respective modes of vibration 
and thus to determine their degeneracies. This has been done (Dayal, 1944 5) 
on the basis of the theoretical formule of Bhagavantam (1943) and of 
Chelam (1943) with the help of three force-constants, two of which, namely 
(K and K"), arise from the mutual displacements of the neighbours and next 
nearest neighbours respectively, while the third Ka is involved in the changes 
of the various valence angles. In the absence of any experimental knowledge 
of the frequencies in the case of silicon, they have been directly calculated 
from a set of three force-constants appropriately assumed. The force 
constants chosen are 


K =1:-2 x 10° dynes/cm. 
K” =0-36 x 10° dynes/cm. 
K, = -- 0:07 x 10° dynes/cm. 


a 


The calculated values of the frequencies with their respective degeneracies 
are given in Table I. These degeneracies refer to the vibrations in a supercell 
which is the unit of the vibration pattern and has twice the linear dimensions 
and eight times the volume of the unit cell of the crystal lattice. Each degree 
of freedom of this supercell has a degeneracy of N/16 in the whole lattice 
arising out of the translational symmetry. 


The choice of the three force-constants, it should be remarked, has not 
been made in an arbitrary fashion, but is based on the experimental data. 
Stitt and Yost (1937) have found that the frequency of vibration of the two 
SiH, groups against each other in disilane (Si,H,) has the value 434-5 cm, 
from which they deduced 1-7 x 10° dynes/cm. as the force—constant for the 
Si—Si bond. The force-constant for the C—C bond in the analogous carbon 
compound, namely ethane, has been found by Bonner (1937) as 4-3 x 10 
dynescm. The force constant K in diamond was found by us (i.c.) to be 
smaller than that in ethane and could be obtained from it on division by 1-4. 
We have, therefore, obtained K in the case of silicon by reducing the Si—Si 
force constant for disilane in the same ratio. The force constant K” has 
been approximately obtained from the observed compressibility. Nagendra 
Nath (1935) has shown that the compressibility 8 of diamond is given by the 
formula 
7 3d 
~ (K + 8K") 
where d is the lattice constant. The compressibility of silicon calculated 
from this formula and the assumed force constants is 0-39 x 10-12 cm.?dyne, 
which compares well with the observed value 0-32 x 10-1? (Richards, 1915). 
The third force constant Ka was found to be 1/16th of the force-constant K 
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in the case Of diamond. The ratio of K to Ka is the same in the present 
case also, but whereas Ka was positive in diamond, it has a negative sign 
here. We have found (Dayal, 1944) that the deformation force-constant 
Sis alwavs negative in the case of metals. A negative sign for Ka is thus 
in keeping with the metallic character of silicon. 


TABLE | 


Vibration Frequencies of Silicon 























Frequency 
No. Mode* hy Degeneracy 
Cm" k 
vy F, 321 460 | 3 
Ve H 1 463 660 | 6 
V3 H, 190 273 6 
Vs Hy 523 750 6 
Vs K; 563 805 4 
Ve K, 363 520 4 
Va My -— . 620 8 
Vg M; 134 192 g 
| } 








* The symbols F2, etc., denote the group theoretical symbols used by Bhagavantam (l.c.). 
3. Evaluation of the Specific Heats 


In our work on diamond we had represented the residual 1/16th of the 
total degrees of freedom, not covered by the eight optical frequencies, by a 
Debye function. In order to take the smaller number of degrees of free- 
dom into account, the usual characteristic temperatures was however, divided 
by 4/16. We have used the same method here. The elastic constants being 
unknown, they were evaluated from the above force constants by the use of 
Nagendra Nath’s formule (/.c.). The values obtained are: 

Cy, = 2°89 x 10% dyne/em.? ¢c,. = 2°34 x 10"? dyne/cm.? 

Cy, = 5:54 x 10" dyne/cm.? 
The characteristic temperature was obtained from the average value of 1/v3 
(v = velocity of elastic waves) along the cube axes and the cube diagonals, 
the statistical weight in taking the average being 3 and 4 along these direc- 
tions respectively. The calculated @ is 530 and the reduced value is 210. 
The formula for C,, then becomes 


c= [dee (i) +207) 


where E and D stand for the Einstein and the Debye functions respectively. 
n, is the number of degrees of freedom of the supercell corresponding to the 
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frequency v, and is given in the last column of Table I. The calculated 
atomic heats are given in Table IT. 


The specific heats C, of this substance have been measured by Nernst 
and Schwers (1914) at a few temperatures between (20-90)° K., by Anderson 
(1930) from (61-296)°K., and by Magnus (1923) from (138-1,000)° K, 
Nernst and Schwers’ values are considerably (2 to 30%) lower than those 
extrapolated by Anderson from his own data. We have shown only the 
latter’s extrapolated values in Table II in this region. (C,—C,) has been 
taken from our earlier work (Dayal, 1941) where it was calculated from the 
approximate formula (C,—C,) =AC*,T. A was obtained from the 
linear coefficient of expansion (7-15 x 10-*), and the compressibility 
(-325 x 10-1*), the values having been taken from the International 
Critical Tables for 23°C. Since, however, the linear coefficient of expansion 
as measured by different workers varies between (2-8 to 7-3) x 10-* (L.¢. 
Tables, Vol. I, page 104), the correction is extremely uncertain. The un- 
certainty in the values of C, due to this cause is only important at very high 
temperatures and may be as high as 3-4% in the neighbourhood of 700°K. 


TABLE II 
Atomic Heats of Silicon in Cals./Deg. Gm. atom 














i 
Temp. Optical Elastic | Total Experimental Literature 
“K. frequencies spectrum | Cy Cy 
20-1 -007 -025 | -032 -042 Anderson 
(extrapolated) 
28-2 -055 -062 | *117 +123 os 
39-7 +224 *125 | +349 -336 ”» 
53°1 +492 -190 | +682 -630 a 
61-2 664 +222 | -886 “811 Anderson 
(experimental) 
65+1 -743 233 | -976 -974 v 
75-0 +949 +259 | 1-208 1-160 ” 
98-8 +435 301 1-736 1-695 % 
-809 +318 | 2°127 2-152 ” 
-067 +324 2-391 2-429 9 
+384 +345 2-729 2-740 a 
*631 +345 2-976 2-993 oa 
+845 +346 | 3-191 3-156 Js 
-102 +349 | 3-451 3 +496 ‘a 








+262 +353 3-615 3-684 99 
*473 +354 3-827 3-841 ” 
+753 -358 | 4-111 4-076 + 
918 +359 4-277 4-335 ” 
+173 +362 ' 4-535 4-566 hy 
591 +366 4-957 4-972 Magnus 
*797 -367 5-164 5-094 | 99 
+128 -370 5 +498 5+298 * 

| 5+672 . 
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In conclusion, the author wishes to record his grateful thanks to Pro- 
fessor Sir C. V. Raman, Kt., F.R.S., N.L., for the interest taken by him in this 
work. 

4, Summary 

The specific heats of silicon have been evaluated on the basis of the 
Raman theory of crystal vibrations. The eight discrete frequencies indicated 
by the theory have been calculated from force constants based on the Raman 
effect data of disilane and the observed compressibility of silicon. The 
residual 1/16th of the total degrees of freedom of the crystal, not covered by 
the eight discrete frequencies, have been assigned to the elastic spectrum and 
their contribution to the specific heats has been evaluated by a Debye 
function with a correspondingly reduced parameter. This has been obtained 
from the elastic data deduced theoretically from the force constants. There 
is an excellent agreement between the calculated and the observed values of 
the specific heats throughout the entire range of temperature (20-1,000)° K. 
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/. Introduction 


DIAMOND has been the subject of numerous studies at this Institute during 
recent years. Nayar (1942), Miss Mani (1944) and R. S. Krishnan (1944) 
have investigated the vibration spectrum of its crystal lattice, and the latter 
has identified 1332, 1284, 1248, 1149, 1088, 1103, 784 and 565 as its charac- 
teristic frequencies expressed in wave-numbers. The presence of so large 
a number of discrete optical frequencies in a crystal with only two non- 
equivalent atoms in the unit cell has been explained by Sir C. V. Raman 
(1943) on the basis of a new theory of the vibrations of crystal lattices, 
Bhagavantam (1943) and Chelam (1943) have obtained expressions for the 
normal modes of the diamond lattice on the basis of this theory, and Dayal 
(1944) has correlated these frequencies with the known specific heat data of 
Pitzer (1938) below 273° A and of Magnus and Hodler (1926) above this 
temperature. Nayar (1941) has studied the scattering of light in diamond 
over a wide range of temperatures and finds that the principal Raman line 
with the frequency shift of 1332-1 cm.-! at 25°C. remains sensibly sharp 
over the whole range of temperature from 85° A to 1130° A, the frequency 
shift at first diminishing slowly and then more rapidly with rising temperature, 
being 1333-8 and 1316-4 respectively at the two extremes of the range. 
These frequency changes are evidently determined by the thermal expansion 
of diamond, and although the data regarding the latter are rather scanty, it 
appears worthwhile to review the existing theories in the light of the observed 
facts. 





2. Theoretical Basis of Griineisen’s Formula 


Griineisen was the first to suggest that the co-efficient of expansion of a 
crystal is related to its specific heat, and showed this to be the case for several 
metals. The relation has come to be regarded as a fundamental one in the 
theory of the solid state, and severai proofs of its validity have been suggested. 
Griineisen himself (1912, 1926 a) assumed a law of force between the atoms 
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in which the attractive and repulsive forces varied as the mth and nth power 
of the distance between the atoms and was thus led to his well-known relation 





i a (1) 
Xe p—mtnts3 Xo ’ 
ae 
_ 1 dV =. Fd .. Xo 
or a, a (2) 
(neglecting the expression — me uth vy E which is very small). 
0 
¥, is here the compressibility at zero temperature, E the energy, 
—dlogv,. . 
i eev a is the thermal expansion coefficient, C, the atomic heat 


and V, the atomic volume. Debye, on the other hand, deduced from 
thermodynamic considerations the equation of state for monoatomic iso- 

d log v 
~ dlog V 
energy. From this equation, Griineisen’s relation follows as a first approxi- 
mation. 


tropic bodies in the form ( p+ *) ¥=- E, E being the thermal 


Dynamically, a crystal can be approximated by a set of 3N vibrations 
with frequencies v,, v2....vs,, there being m atoms in the crystal. If the 
anharmonicity is small, it can be shown (Slater, a) that the coefficient of 





. ° . _ kek ‘ ae ae _ hy; 
thermal expansion is given by a = V, Ly; ( —1)3 where x = kT and 
y= = ve Griineisen assumed that all the »;s were equal to each 


other in a regular monoatomic lattice, so that y may be regarded as an 
empirical constant. We have then Griineisen’s relation in the form already 
stated. 


The justification for Griineisen’s formula (Slater, 5) comes from the 
9 log Vinx 
3 log V ’ 
spectrum is entirely determined by the limiting frequency v,,,, and if this 
frequency changes, all other oscillations change their frequences in the same 
ratios. Griineisen (1926, c) found support for his formula in the case of 
diamond. Taking the Debye temperature @ to be 1860 as found by Nernst 
and Lindemann from the specific heat data, he found an agreement between 
the observed and calculated values of the thermal expansion. We shall, 
however, now re-examine the situation in the light of Nayar’s spectroscopic 
data, 


Debye theory giving y = — according to which the frequency 
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3. The Failure of Grijneisen’s Formula 


d log v 1 Av 1 
We have y= — > log V “~> i-th. 
efficient of thermal expansion and Av the diminution in frequency » in 
the termperature range T, to Ts. 
“ Av=ya(T, —T,)» (3) 
We use equation (2) to find a the mean coefficients of thermal expansion 
and equation (3) to find the diminution in frequency of the line 1332-1 (25°) 
in diamond with the rise of temperature. We need the mean values of the 
specific heats at constant volume. The latest specific heat data (C,) are by 
Pitzer below 273° A and by Magnus and Hodler above this temperature. 
These have been plotted on a curve and from it C, was obtained by finding 


where @ is the mean co- 


Ta 
the area Fs C,, dT and dividing it by (T,— T,). The bulk modulus of diamond 


T1 
at 20°C. has been determined by Adams (1921) and Williamson (1922). 
Their values are respectively 6-25 x 102 and 5-56 x 10!2. Griineisen 
found a value of y = 1-1 by using the first value. We shall also use this 
value of y in the present calculation. The observed values of the mean 
coefficient of thermal expansion are those of A. Dembowska (quoted by 
Griineisen, c). The results of the calculation are shown in Table I. We 























TABLE I 
one Mean value of thermal expansion a | Av 
v . 
Temp. Range A Mean value Ny 

of C, Calculated x10, | Observed x 10, 7 ynantery 
84+8- 194-1 +2333 -502 | +54 7 é 
194-1-— 273-2 *866 1-860 1-74 ie one 
273 +2— 296-2 1-351 | 2-907 2-91 mA iia 
296+2- 328 1-605 3-452 3-51 ais Ss 
328 —- 351-1 1-870 4-02 4-35 Si a 
300 - 400 1-963 4-224 Aes *62 1-2 
400 —- 500 2-921 | 6-284 92 | 164 
500 — 600 3-691 7-939 1-16 16 
600 —- 700 4-166 8-962 1-31 1-9 
700 - 800 4-466 9-607 1-41 2-0 
800 —- 900 4-710 10-14 1-48 2-2 
900-1000 4.913 10-57 1.55 5.4 

13D) | 255 


| 


1000 -1100 | 5-041 10-84 a | 
| | | 








use Xag3 instead of X, as the difference would be very small for a crystal of 

extremely low compressibility like diamond. For higher temperatures 
+n+3 Xp 

6 ’B, E, 


Griineisen includes in equation (1) the correction term 
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although he neglects it for the region in which the thermal expansion data 
are available for diamond. 
According to Griineisen, = should lie between 2 and 3. 


ve = 457K. calories for diamond and Ej,y,9=—3-3K. cal. Using the 
0 

maximum value 3 for = — 
in the calculated value of a in the highest temperature region. Further, the 
mean values of c, should have been used in the above calculation. These 
however differ from the C, values by only 1% at the highest temperatures. 
The calculated value at high temperatures of Av is thus only very slightly 
(to the extent of 3%) affected by using C, and neglecting the correction 
term. We find from Table I that though the Griineisen constant 1-1 gives 
a fairly satisfactory agreement between the observed and calculated values 
of a, it completely fails to fit Nayar’s spectroscopic observations. 


we get a maximum increase of 2% 


4. Thermal Expansion on the New Lattice Theory 


The lattice spectrum of a crystal according to Sir C. V. Raman exhibits, 
in general, (24p — 3) monochromatic modes of vibration, p being the 
number of non-equivalent atoms in the basis cell. These are the modes of 
vibration of the atoms in a super-cel! having the twice dimensions and 
eight times the volume of the unit cell. The three missing degrees of 
freedom represent the elastic spectrum of the solid. Accordingly, the mean 
coefficient of thermal expansion may be written as 
_ Xo 3RJ “Ee d log v; 


_\ 0 log » 
“aut.” Coy log V + 3 Dm) d log vy 


In this equation C,, represents the contributions of the monochromatic fre- 
quencies v; to the specific heat C, and o; is the statistical weight of frequency 


a 


v; and —- i y; is its volume dependence. It has been shown thermo- 


dynamically that if we assume y to be the same for all frequencies, y is inde- 
pendent of temperature even at low temperatures (Mott and Jones). We 
take therefore y,; to be independent of temperature. For obtaining the mean 
values of ermal expansion, we have to form energy differences so that 
_ _ Px(U) — Pr(U) 
Ts | T; 
_Xq 3RJ 


where PU) =} a ["2 ” py MEU ++ 3 ee Dw) 


E (U) and D(U) being the Einstein and Debye energy functions. 
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The frequencies of the diamond lattice are: 1331 1284 1248 1149 1088 
3 8 6 6 4 
1013 784 565, the number below each frequency indicating its statistical weight 
4 6 8 

according to Dayal. The Debye limiting frequency has been calculated as 
vy = 573 and has weight 3. We have no knowledge of the volume dependence 
of these frequencies, and we have therefore to make some reasonable assump- 
tions regarding their y values which will fit Nayar’s data as also the observed 
values of expansion cofficient. The higher frequencies 1332, 1284, 1248, 
1149, 1088 and 1013 represent modes of vibration in which the distances 
between neighbouring carbon atoms alter and are determined by the binding 
forces between them, while the lowest two frequencies 784 and 565 do not 
involve any variation of the distance between the carbon atoms linked with 
each other. We divide the frequencies into these two groups and assume 
that all the frequencies in each group have the same value of y. Further, 
there is no reason to suppose that the transverse modes of vibration in the 
elastic spectrum would have the same volume dependence as the longitudinal 
modes, since the nature of the two modes is essentially different. We there- 
fore divide the elastic spectrum into two branches, viz., the longitudinal and 
transverse branches, the former with statistical weight one and the latter with 
weight two. We find by calculation v, = 842-5 and v, = 514. These values 
give good agreement with the observed specific heat data. We include », in the 
group of higher frequencies and assume the same value of y for it and », in the 
group of lower frequencies with a different y. We thus need only two different 
y values for the crystal, y, for the group of higher frequencies and y, for the 
group of lower frequencies. We accordingly write P(U) in the form 


P(U) = ens [ye {3E, (U) + 8E, (U) + 6E;(U) + 6E,(U) + 45, (U) 


~ Vy 48 
+4 E,(U) + Dg (U)}+ y {6E,(U) + 8E,(U) + 2D, (U)}] 


where E, (U), etc., are Einstein energy functions for the 8 monochromatic 
frequencies and D,(U) and D,(U) are the Debye energy functions for the 
longitudinal and transverse Debye frequencies. 


We have chosen y, = 1-402 and y, = 0-834 to fit Nayar’s data as also 
the thermal! expansion coefficients. In choosing two constants which fit 
two separate sets of data, we are not really assuming two arbitrary constants 
which can be varied at will. Moreover. it is prima facie evident that the 
modes of vibration which involve change in the C — C bond length would 
contribute to the thermal expansion in a larger measure than those which 
do not. 
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The Thermal Expansion of Diamond 


As diamond is a homopolar crystal, we may use the method given 
by Slater (c) for finding y. 


3 

y= So} an ot = Ne 
where L is the heat of dissociation, c = 8/3*, ra the C—C distance, x» the 
compressibility, and N the Avagadro number. For diamond; L = 72-6 
K. cal. per gram atom (Fajans, 1920, 1922) and rp = 1-54 x 10-8 which gives 
= 1-17 if yg =0-180 x 10-7%. It will be noticed that this is about the 
same as the value of y assumed by Griineisen and lies midway between the 

values of y. and y, found above. 


5. Calculation of Thermal Expansion Co-efficients 


We use y, = 0-18 x 10-!2 which is the value obtained from the latest 
work of Williamson. The results of the calculation are shown in Table II. 











TABLE II 
| a x 10° | Av 
Temp. Region ! ! 

Calc. Obs. | Cale. | Obs. 

84-8- 194-1 “49 +54 | = we 
194-1- 273-2 1+85 1-74 a et 
273-2- 296-2 | 3-15 2-91 ee Re 
196-2- 328 3-64 3-51 ss * 
oe. | 4-31 4-35 ie - 
194-1- 296-2 2+19 - “41 8 
296+2- 400 4-52 -88 1-2 
400 ~ 500 6-82 | | 1-28 1-35 
— 600 8-72 | 1-63 1-60 
600 - 700 10-16 | | 1-90 1-85 
700 - 800 | lies | 2-08 2-0 
800 - 900 12-13 2+24 2-2 
900 -1000 12-77 2+36 2-4 
1000 -1100 13-23 3+20 3+3 











We find from Table II that the thermal expansion data and the spectro- 
scopic data of Nayar are both reasonably well fitted. There is some diver- 
gence between the observed and calculated values of A v in the low temperature 
region, but this is most likely due to the difficulty of measuring the very 
small frequency changes in this region. In the high temperature region 
the agreement is remarkably good. 


Using the same values of y, and r,, we have also calculated the co- 
efficients of thermal expansion at 20°C., 40°C., 50°C. for which the data 


are available from the J. C. Tables and Physico-Chemical Tables, Vol. I, 
p. 17, 








Bishambhar Dayal Saksena 











Temp. °C. a x 10%(Calc.) aX 10* (Obs.) | Observer 

{ 
r 

20 3-02 2-7 | < < 5 

2-91 | H.B.P. 

40 3-64 3-54 Fizeau P.C.T. 

3-88 3-96 | Fizeau P.C.T. 
} 








(Griineisen gives the values of a at 20°C. as 2-91 x 10-*). The calculated 
values show good agreement with the observed values of Fizeau. 


It is to be noted that the values of the co-efficient of thermal expansion 
calculated by us for the region of higher temperatures (vide Table II) are 
greater than those calculated on the basis of Griineisen’s formula (Table 1). 
This is a point which can be tested out experimentally. 


In conclusion the author thanks Professor Sir C. V. Raman, for his 
kind suggestions and keen interest in the work. 


6. Summary 


It is shown that the well-known Griineisen’s law of thermal expansion 
does not fit in with Nayar’s observations of the change of the Raman fre- 
quency shift in diamond with temperature, although it gives a fairly good 
agreement between the observed and calculated values of thermal expansion 
in the lower ranges of temperature for which alone data are at present 
available. The assumption of a single Griineisen constant for diamond is 
not in agreement with the facts and does not appear to be justified. The 
expression for thermal expansion on the basis of the new lattice theory of 
Raman has been written down. The frequencies have been divided into two 
groups, one of higher and the other of lower frequency and the volume 
dependence has been assumed to be the same for all the vibrations in each 
group. By assuming that y is different for these two groups, but the same 
within each group, it is possible to fit both Nayar’s spectroscopic data and also 
the thermal expansion co-efficients in the range for which there are observed 
data. The values of the co-efficient of thermal expansion calculated from the 
new expression are definitely larger at high temperatures than those to be 
expected on the basis of Griieisen’s formula. The assumption of two differ- 
ent ys for the two groups of vibrations also appears prima-facie justifiable, 
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7. Introduction 


DaRWIN (1914) and Prins (1930) have investigated the intensity of reflection 
of X-rays by crystal planes, and have derived formule for the variation of 
intensity with the angle of incidence. The special case of a non-absorbing 
crystal was studied by Darwin, who found that the reflection is total over a 
small range of angles of incidence, outside which it falls away rapidly. 
Prins has extended Darwin’s formule to an absorbing crystal. In essence 
the method is identical in both cases; it depends on the formation of a set 
of difference equations, and the solution of these for the case when the number 
of planes in the crystal is infinite. An excellent resume of the Darwin-Prins 
analysis is given in the book, “‘ X-rays in Theory and Experiment,” pp. 365 
to 391, by Compton and Allison (hereafter referred to as C and A), whose 
notation we shall follow in this paper, on account of their compact nature. 


An alternative method would be to solve the equations for the case when 
the number of planes is finite and then to extend the results for the infinite 
crystal. This procedure is more instructive in that it shows the manner in 
which the principal maximum is built up. In fact. it is found, as will be 
shown later, that, outside the region of perfect reflection indicated by 
Darwin’s theory, the curve is not smooth as shown by his, but that there are 
an infinite number of maxima and minima, which, however, are so close 
together that it will not be possible to detect them. This important feature 
is altogether missed in Darwin’s theory. 


The author has actually considered such a finite stratification in an 
earlier paper (1942). There, a most general case of a periodic stratification 
was considered, in which the variation in the properties of the medium was 
supposed to occur periodically, the nature of this variation being unspecified. 
The theory can easily be applied to crystal planes, and then it can be extended 
to the case of a crystal of infinite depth. This is done in this paper, and the 
results obtained are substantially the same as those of Darwin and Prins. 
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It is however found that Darwin’s theory leads to a different formula from 
that in the present paper for the region outside that of perfect reflection, 
which obviously arises from the failure of his theory to take cognisance of 
the fluctuations of intensity which occur in this region. Incidentally, it is 
found that a certain assumption made by Darwin and Prins is not altogether 
justified. The present method is free from any such, and is consequently 
more rigorous. 


2. General Formule 


Suppose that we have a medium consisting of 7 stratifications, and that 
a beam of electromagnetic radiation is incident on it. Let T, be the ampli- 
tude of the incident beam at the beginning of the first stratification. It is 
evident that there will be a large number of multiple reflections, and conse- 
quently there will be two streams of energy in the medium, one building up 
the transmitted beam, and the other the reflected beam. Let T,, R,; T;, Ro; 
..3 Tys, represent the electric vectors in the two streams at the beginning 
of the first, second, etc., stratification. If now we denote by r and ¢ the 
complex reflection and transmission coefficients of a single stratification, 
on which the radiation is incident at the angle concerned, then the following 
difference equations are obtained: 


R, =7T, +tRwa; T, =tT,4 +R, (1) 
Eliminating the R’s and the T’s successively from these, 
Tha = yT, ai Tyan R,4 bss YR, Rys (2) 


where y = (1 + ¢2 — r*)/t. These can be solved by making use of the 
condition that R,,,, =0. One obtains in this way 


RV Tra =C/DAO); Ti/Taa =f. O)/t —-faly) (3) 
where Liar = OD = yea Om oN ao (4) 


The series appearing as f,,(y) can be summed up by putting y = 2 cosh 8, 
and comes out as =sinh n§/sinh 8 (Jolley, 1925). Also, making the 
substitution 
ee ee 
sinh 8 sinha K 
where K can be evaluated to be equal to sinh (a-+ 8) from the relation 
cosh 8B = (1 +f? — r?)/2t, one finally obtains, 


Ry Tres T, 


sinh nB sinha sinh (a+ np) (5) 
where cosh 8 = (1 + £2 — r?)/2t and cosh a =(1 — ¢? + r?)/2r (6) 
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From (5), the reflecting power, R, of the stratified medium is given by 
2 
R= Ri | =! e0sh « + sinh « coth nB |? (7) 
This is our fundamental equation, from which the intensity of the reflected 


beam can be calculated by substituting the values of a and 8. It may be 
noted that, in general, both a and f are complex. 


4 


3. Derivation of the Formule for an Absorbing Crystal 


In order to determine the values of a and 8, we must know those of 
r and t. For these, we make use of the relations, first derived by Darwin 
and reproduced in C and A, pp. 365 to 376. If we take the boundary 
between adjacent stratifications as the planes midway between the atomic 
planes, it is easy to show from the above formule in C and A that 


r = —ise- 2!2 and t =(1 — io)e* 82 (8) 


where s and o are two constants depending on the nature of the planes, 
k =2mn/A, and 5 = 2d sin y, d being the distance between the planes, ¥ the 
angle of incidence, and A the wavelength of the X-rays. The quantities 
o and s are given by 

(A sin ¢/2nd)o =p, — 1 = — (8 +8) (9) 
where », is the complex refractive index of the crystal, whose real and 
imaginary parts are, respectively, (1 — 5) and 8, and 

s =f (2p, k) o/f (0, k), 

so that (A sin %/2nd)s = —(8 + i). f (2%, i/f(0.k) = —(A +75) (10) 
where f (2%, k) and f(0, k) are quantities proportional to the atomic struc- 


ture factors for angles 2 and zero. (C and A, p. 376). In the above 
expressions, o and s are to be taken as small. 


Now, we are interested in the phenomena that occur near about the 
Bragg angle, for which k 5/2 ~ mz. Putting k 5/2 = m7 + &, where & is 
small, exp (— ik 8/2) =(— 1)” exp (— i€), and 

r =(— 1)"*" is et, and t =(— 1)” (1 — io)e~t (11) 
From these, treating 0, s and é as small quantities, one obtains 
cosh a = (— 1)”+! (0 + €)/s; cosh B =(— 1)” {2 +5? — (o + €)*}/2 (12) 


In the general case of an absorbing crystal, both o and s are complex 
quantities, so that both cosh a and cosh 8, and hence a and 8, are complex. 
Now, the reflecting power of a finite stratification is given by (7). In order 
to determine the same for an infinite one, we must put n —>co in the right 
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hand side of (7). If B is complex, it is easily shown that Lt coth n8 = 1, so 
that for an infinite stratification, we 


Lt kA = cosh a + sinh a (13) 
n->co 1 

Now, from (12), cosh a =(—1)”*! (¢ + €)s, so that sinh a = 

tV(o + &)? — s*/s, and 
Lt B= —si(— I" @ +H 2versr—s} (14) 

n>co *1 

This is identical with the Eq. (6-42) in C and A, p. 380, except for the factor 
(- 1)” before (o + €) in our above expression. The difference is easily 
explained, for the amplitudes are measured at points on the planes in the 
derivation in C and A, while in ours, the reference point is at a distance d/2 
above the first plane. It is obvious that this would produce a phase change 
of exp (— ik8/2), which, for small values of &, reduces to (-- 1)” exactly the 
factor by which the two expressions differ. It is thus seen that Prins’s formula 
can be derived easily from our general formule by substituting the appro- 

priate values of r and f. 


4. The Case of a Non-Absorbing Crystal 


We now consider the case of a non-absorbing crystal. Here, 8 and bd in 
Eqns. (13) and (14) are equal to zero, so that both o and s are real. Hence, 
both cosh a and cosh 8 are real, and two separate cases have to be considered, 

* 


viz., 
case (i) | cosh a | < |, for which | cosh 8 | > |, and 
case (ii) | cosh a | > |, and | cosh 8 | < |. 
From (12), |cosha| =1, if(o+@)= +s,0r§=—os. 


Hence, if —o —s <& < —wa +5, case (i) is operative, so that sinh? a is 
—ve and sinh?8 is-+-ve. If one writes a =a, +ia,, and B = 8, + i, then 


ms. sii h bar sin?a, 

R =| COS a2 + / sin a, cot np, |? =1 + sinh? 78, (15) 
Therefore, as n is increased, the intensity corresponding to any value of ¢ 
within this range steadily increases. When n—co, sinh? nB,—oo, and 
R—1, ie., the reflection is total throughout the whole range of values of & 
from —o—s to —o +8 for an infinite stratification. This region may be 
called the principal maximum, and its width is given by A & = 2s. 

A3a 
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Outside this range, case (ii) is operative, so that sinh*a is +ve and sinh?¢ 
— ve. Hence, in this case, 


1/R = 1 + sinh®s,/sin? 8, (16) 
In this region, a number of maxima and minima occur corresponding to 
nB, = sx + 1/2 and sz respectively. These may be referred to as subsidiary 
maxima, and they all lie evidently on the curve R =sech*a,. When the 
number 7 increases, the intensity corresponding to any value of & in this 
region does not tend to a limit, but only undergoes fluctuations. Also, the 
subsidiary maxima approach closer and closer together. When n actually 
tends to infinity, the subsidiary maxima would be so close together, and the 
fluctuations in the intensity would be so rapid, that they would not be detected. 
The observed intensity would only be an average over a cycle of changes, 
and the observed value of R would be given by 


| O+® v7 in n8 
Ra 2 f RAM (scene, esinnta, AMA (17) 
6 U] 


sin? nB, + sinh? a, 


The integration can be carried out, and leads to the remarkably simple result 
that 
R,,, = 1 — | tanha, | (18) 

Now, the boundary between the two regions is obviously given by sinh a, = 0, 
for which tanh a, = 0 and R,,, = 1. Thus, in the second region, R,,, 
starts from unity and quickly drops down as the value of & is taken farther 
and farther away from a principal maximum. These features are all very 
similar to those predicted by Darwin’s theory. In fact, the reflection curve 
is not symmetric about ¢ =0, but € = —a, or, k 8/2 =mn —a. The 
value of this shift caused by refraction, as also of the width of the region 
of perfect reflection can be obiained in angular dimensions. If 6, be the 
Bragg angle defined by k 5 2 = kd sin 6) = mr, and if @ be the angle of 
incidence corresponding to the centre of the Jiffraction pattern, then 

0 — 0g = — o kd cos Oy = 5 sec Oy cosec 44 (19) 
Also, if the angular width of perfect reflection be Au, then 

Au =s/kd cos 9 = A sec % cosec 4 (20) 
Both these are identical with the relations derived in C and A, pp. 388 and 
389. . 


5. A Note on Darwin’s Method of Solving the Difference Equations 


In the preceding section, it was shown that the reflection coefficient in 
the region outside the principal maximum is R = | — |tanha,|. Darwin’s 
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formula, put in terms of our symbols, gives R = e-!?™! in the same region, 
The disagreement needs an explanation. It arises from a certain assumption 
made by Darwin in solving the difference equations. In terms of our symbols 
Darwin puts T,’/T, = x, and obtains R,/T, in the form 


R,/T, =(r? — ¢2 + tx)/r (21) 
The problem thus reduces to the evaluation of x, and for this, Darwin makes 
the assumption that T,,, =T, By substituting this in (2), he gets 
x =e*f=cosh B+ sinh 8 and R,/T, =cosha + sinha. (22) 

For the absorbing crystal, the same procedure is followed (C and A, p. 378). 


It is to be noted that such an a priori assumption cannot be justified, 
since we are actually trying to find the values of the quantities T, and R.. 
A more serious objection to this assumption lies in the fact that it leads to an 
internal contradiction in the theory. Taking the case where it disagrees 


with the present theory, the value of x here is exp (+ i8,), so that |x| =1. 
Hence, | T,,, | =| T; |, so that Lt | T,| =| T, |. Thus the transmitted in- 
8—->co 


tensity is equal to the incident intensity and the reflection must be zero. How- 
ever, by substituting the value of x in (21), Darwin actually gets R =| cosh a 
+ sinh a |? = exp (+ 2a,). Thus proceeding in two different ways from the 
same assumption, one gets different results, so that the assumption itself 
cannot be correct. 


In fact, by making use of the method of the present paper, it can be 
proved that the assumption holds asymptotically for the first few planes 
(i.e., for small values of s) in Prins’s case, and in the region of perfect reflection 
in Darwin’s case. In other words, the value of x as defined by the equation 
x =T,/T, is actually equal to exp(+ 8) in these cases. However, the 
assumption is not generally true, and is invalid in the region outside the 
principal maximum for a non-absorbing crystal. Details are not given for 
want of space. 


In conclusion, I wish to express my thanks to Prof. Sir C. V. Raman for 
the many helpful discussions I had with him during the investigation. 


Summary 


In this paper is described a new derivation of the formulz for the reflec- 
tion of X-rays by perfect crystals, which forms an alternative approach to 
the problem to that adopted by Darwin and Prins. It consists in obtaining 
a solution of the difference equations which occur in the problem for a 
crystal containing a finite number (n) of laminations. The case of a crystal 
Aa 
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of infinite depth is obtained from this by proceeding to the limit when n—>.o. 
The formule thus obtained are identical with those of Prins for an absorbing 
crystal, while for a non-absorbing crystal, the width of the region of perfect 
reflection is in agreement with Darwin’s value. However, there is gq 
difference in the formule for the variation of intensity with angle outside 
this region. This has been shown to be due to an assumption made by 
Darwin which is not altogether justifiable. 
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AN earlier paper! described attempts to synthesise 9-8-aminoethyl-9: 10- 
dihydro-phenanthridine and a series of possible sympathomimetics derived 
from some carbocyclic ring systems. The latter had the repeating pressor 
active, atom-patterns, Ar.C.C.N.—as features of their molecular archi- 
tectures and, on the basis of von Braun’s conception of the mode of action of 
the intensely sympathomimetic methyl amino hydrindene*, were entitled 
to consideration as f-phenylethylamines many times over. The present 
communication is concerned with attempts to synthesise further 9-8-amino 
ethyl-9: 10-dihydro-phenanthridine and N,N’-trimethylane bis-tyramine 
which also comes under the same category. 


The scheme for the synthesis of the 8-aminoethyl phenanthridine deri- 
vative could not be pursued owing to failure at the cyclisation stage of 
cyanacetamido diphenyl to 9-cyanomethyl phenathridine, a necessary 
intermediate. However, success was encountered in the scheme for the 
synthesis of the trimethylene dis-tyramine. 4-Methoxy phenyl ethylamine*® 
was transformed in good yield to its N-benzenesulphonyl derivative‘; the 
potassium salt of the latter gave with trimethylene bromide the correspond- 
ing alkali insoluble di-amide, which on hydrolysis, furnished the desired 
N,N!-trimethylene bis-tyramine. 


OCH; 
(CH,)3Br OCH; 
| | > CoA foes 
ys CH *CHg*N -(CHg)3*N+CH2:CHo+CgHy 


l 
ha oe 
CH»+CHy+NH -SOg'CgHs SOs°CoHy SO2 -CoHs 


— a os 
H0o-C \-cn, *CHg*NH°CH2:CHg ‘CHg: NHCH,-CH2—-7 


\ 
K P >=-OH 
aad N, N’—Trimethylene dis —tyramine * 


N,N’-Trimethylene bis-tyamine has recently been subjected by Capt. 
K. Venkatachalam of the Madras Medical College to preliminary biologi- 
cal tests in respect of its pressor activity and compared to tyramine in the 
spinal cat, The tests, so far conducted, reveal that while its intensity of action 
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is nearly equal to that of tyramine, trimethyle ie bis-tyramine is superior to 
the standard sympathomimetic in the matter of its duration of action. 


Experimental 

2-Cyanoacetamido diphenyl—A mixture of 2-amino diphenyl (4g), 
and ethyl cyanoacetate (2-8 g.) was heated under reflux at 200-10° for 2 
hours and cooled; the solid (2-5g.) which separated was washed with 
dilute hydrochloric acid followed by water and crystallised from water, 
It was obtained as colourless needles, m.p. 126° (Found: N, 11-3. C\;H,,ON, 
requires N, 11-9 per cent.) 

N,N’-Trimethylene bis-tyramine hydrochloride. The postasium salt 
of N-benzenesulphonyl tyramine-methyl ether was prepared in situ by add- 
ing the calculated quantity of potassium to a hot solution of the benzene 
sulphonamide (7g.) in xylene (15c.c.) Trimethylene bromide (2-4g) 
was then added and the mixture refluxed for 4 hours. The reaction mixture, 
after pouring into water and steam distillation to remove xylene completely, 
furnished the trimethylene diamide. The crude product was freed from 
the unreacted starting amide by washing successively with dilute alkali and 
water. The purified diamide was then heated with concentrated hydro- 
chloric acid (10 c.c.) in a sealed tube at 160-170° for 6 hours so as to obtain 
hydrochloride of N,N’-trimethylene bis-tyramine. 


The crude, dark coloured product of the pressure reaction was dissolved 
in absolute alcohol and fractionally precipitated with anhydrous ether, the 
first oily portions being rejected. The partly purified hydrochloride (1-4 g.) 
then recrystallised from alcohol-ether in colourless needles, m.p. 149°, 
(Found: N, 7-8; Cl, 18-6. C,9>Hss0,Ne, 2 HCl requires N, 7:2; Cl, 19:0 
per cent.) 

I am thankful to Professor P. C. Guha, D.Sc., F.N.I., for his kind 
interest in this investigation and to the Government of Madras for a Reserch 
Scholarship. 

Summary 

In continuation of an earlier line of work, synthesis has been effected 
of N,N’-trimethylene 4is-tyramine which fulfils the requisite structural 
specifications for sympathomimecity. The results of subjecting this to 
preliminary biological tests in respect of its pressor action in the’ spinal cat 
are also indicated. 
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